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We show that the neutron scattering spectra observed in CaFe2As2 by Zhao et al 
(arXiv:0903.2686vl) and the highly unusual spatially anisotropic exchange constants in the iron 
pnictides derived by Han et al (Phys. Rev. Lett. 102, 107003 (2009)), using electronic structure 
calculations, can be explained by assuming a role for orbital order in these materials. We write 
down a simple model Hamiltonian with tertagonal symmetry, whose spin- wave spectra describes the 
observed dispersion relations. We further argue that these materials have orbital selective Mott tran- 
sition, which is driven by superexchange between neighboring iron atoms. We suggest that reduced 
spin and quasi-one dimensionality and not frustration are responsible for the reduced moments in 
these materials. 
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The discovery of superconductivity at relatively high 
temperatures in the iron pnictide family of materials has 
invigorated condensed matter researchfi In order to de- 
velop a fundamental understanding of high temperature 
superconductivity, it is important to understand the rela- 
tionship and differences between these pnictides and the 
even higher temperature superconductors, the cuprates. 
Both family of materials arc layered quasi-two dimen- 
sional electronic systems, which at stochiometry develop 
antifcrromagnctic order. Upon electronic doping of these 
layers, they become superconducting. One of the key 
differences between the two is that the cuprates, at sto- 
chiometry, are in many respects classic Mott insulators. 
In local density approximation (LDA) based calculations 
they are found to be metallic and only strong correla- 
tion effects drive them to insulating and antifcrromag- 
nctic behavior. In contrast, the parent state of iron pnic- 
tide materials is a metal, whose conductivity increases 
upon the onset of antiferromagnetism^2, This metallic 
antiferromagnetic state is also reproduced in the LDA 
calculations M^^^^^i^ 

Furthermore, the cuprate ions have a single d-hole, 
whereas the iron ions in the pnictides have multiple d 
electrons. Thus the former have spin-half, making up 
a system that has maximum quantum fluctuations and 
possibly exotic incipient quantum spin-liquid behavior 
whereas the local moment in the latter could be much 
larger and thus the magnetism more conventional. Any 
role of strong correlations in the pnictides has also been 
questioned, and this would have a strong bearing on any 
universal understanding of the mechanism of supercon- 
ductivity in these family of materials. 

We focus here on the parent insulating mate- 
rial. Magnetism in these materials has been highly 
dcbated^iSiiiii^ii^iiiii^ii^ Some argue that these mate- 
rials are spin-density-waves driven by nesting of fermi- 
surfaces and they should be viewd as weakly correlated 
itinerant magnets. Others have argued that despite the 
metallic behavior, strong correlations are present and the 
correct starting point for understanding the magnetism 
in these materials is a Heiscnberg type model. Indeed, 



neutron scattering spectra can we be well understood by 
the latter approach. 

The conventional strong correlation view has been that 
the magnetism is highly frustrated by having nearest 
neighbor antiferromagnetic exchange Ji compete with 
second neighbor antiferromagnetic exchange J2 i^iiiii^ii^ 
It is well known that this leads to a (tt, 0) ordered state 
at large J2 as observed in experiments. In this case, 
there is an additional Ising degree of freedom associ- 
ated with the direction of antiferromagnetic and ferro- 
magnetic neighbors, which leads to a finite tempera- 
ture phase transitionj ^^i^^ At, or below this transition, 
small anisotropics and or weak three-dimensional cou- 
pling should lead to antiferromagnetic long-range order. 
The Ising transition should couple to lattice distortions 
and lead to a loss of tetragonal symmetry. Indeed, such 
structural and magnetic transitions are observed either at 
the same temperature, or with the structural transition 
being slightly higher in temperature, as expected from 
theory. The proximity to a quantum critical point near 
J2 ~ Ji/2 can lead to a highly reduced moment which 
fluctuates with small changes in material parametersJ^ 

The key to distinguishing different scenarios is Neu- 
tron Scattering and the measurement of spin-wave dis- 
persion in the full Brillouin zonej ^°i^^ In this respect, the 
Ji — J2 model scenario is also of deep theoretical signifi- 
cance from the point of view of quantum magnetism, be- 
cause the (tt, 0) order is selected by an order by disorder 
phenomenal^ The accidental degeneracy of the classi- 
cal model leads to zero energy states near 4-points of 
the Brillouin zone (0,0), (tt, 0), (0,7r), (7r,7r). Quan- 
tum fluctuations, which lift the accedental degeneracy, 
cause a gap at (0,7r) and (tTjTt), leaving gapless excita- 
tions only at (0, 0) and (tt, 0) as required by Goldstonc's 
theoremji2i2^ 

Recently Zhao et al have performed^i neutron scatter- 
ing experiments on the material CaFe2As2. The experi- 
ments show sharply defined spin waves throught the zone, 
not just at very low energies. The absence of a Stoner 
decay, leads them to strongly favor a local moment type 
picture for the system. However, a big surprise is the 
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finding that the spin-wave energy along the K direction 
is a maximum at (tt, tt). We will show below that this re- 
quires that the nearest neighbor exchange be strong and 
antiferromagnetic in one direction and weak and ferro- 
magnetic in the other. Quantum fluctuations can lead to 
a gap at (tt, tt) but they cannot turn it into a maximum. 
The materials clearly lack tetragonal symmetry and so 
such a result can be considered as resulting from lack of 
tetragonal symmetry. However, one needs to understand 
why a small distortion can lead to such a dramatic differ- 
ence in exchange constants. Indeed, this is exactly what 
the elctronic structure calculations of Han et al using lo- 
cal spin-density approximation (LSDA) find, within the 
tetragonal phase.^^^ These results have not yet found a 
simple explanation. The main purpose of this paper is to 
give a simple scenario by which this happens. 

Rather than trying to write down a complete Hamil- 
tonian for the low energy behavior of this system, we 
proceed in a step by step manner to uncover the physics 
behind the unusual neutron scattering spectra and the 
unusual LSDA results. We begin with a two orbital 
model (which we take as d^^ a-nd dy^) on the square- 
lattice in the XY plane, and assume that the total occu- 
pancy of these orbitals is one. Clues to the origin of such 
a model is discussed later. When there are occupied dxz 
orbitals on neighboring sites along the x-axis, there is a 
strong antiferromagnetic super-exchange Ji between the 
orbitals. Similarly, if there are dyz orbitals occupied on 
neighboring sites along the y-axis, there is strong antifer- 
romagnetic exchange Ji . One can visualize that this is a 
direct exchange between the d-orbitals and is only non- 
zero when the orbitals are oriented towards each other. 
In addition, we assume there is a ferromagnetic nearest 
neighbor exchange Jf and an antiferromagnetic second 
neighbor exchange J2 . One can imagine that both Jf and 
J2 are mediated by the As atoms and are much weaker 
than Ji. The precise variation of J2 and J/ with orbital 
occupation is not important in our model. In general, 
J2 and magnetic frustration does not play a significant 
role in our considerations. So, the full Hamiltonian, with 
tertagonal symmetry, can be written as: 

(1) 

Here Pf (Pf) are one if the dxz {dxy) orbital is occupied 
at i and zero otherwise and Si are the spin operators at 
site i. The sum over < i, j > runs over nearest neighbors, 
and the sum over < i,k > runs over the second neighbors 
on the square lattice. 

The ground state of the model is orbital ordered with 
Pf = 1 and = for all i or vice versa. The 
spin system in this orbital ordered state is unfrustrated 
and has a nearest neighbor antiferromagnetic coupling 



Jia = Ji ~ Jf along X and a ferromagnetic coupling 
Jib = — Jf along y. Thus the system will have Neel order. 
Linear spin- wave theory should be a good approximation 
for the spectra of this model. Using Holstein-Primakoff 
transformations, we obtain the spin-wave dispersion; 

uJk = ^SJ2,J{Al-Bl) (2) 

with 

Ak = I + a ~ f] + l3cosky, (3) 

and, 

Bk = cos(fc:i-)(cos(fc.y) -t- a). (4) 

Here, a = Jia/(2J2), and /3 = J2a/(2J2)- The condi- 
tion for (tt, tt) point to be a maximum along ky is that 
a > 1 and /? be negative and small. A sketch of the re- 
sulting spectra is shown in Fig. 1. Parametrically this 
is the model obtained in the electronic structure cal- 
culations of Han et ali^ And, this is exactly the form 
of the model that gives an excellent fit to the neutron 
scattering spectra, where Zhao et al quote the values for 
CaFe2As2 {SJia = 49.9±9.9meV, SJu = -5.7±4.5meV, 
SJ2 = 18.9±3.4meV),^ 

Because this model has a local Ising like degree of free- 
dom associated with orbital occupation, it should have 
a finite temperature phase transition. However, the en- 
ergetics of the Ising order are mediated by ordering the 
spins. So, it is not clear that the universality class of 
the transition must necessarily be that of the 2D Ising 
model. There could be long-range effective interactions 
and that may lead to a first order phase transition. In fu- 
ture work, such details can be investigated numerically. 
Such an orbital order will clearly couple to lattice dis- 
tortions and a lack of tetragonal symmetry. Note, that 
unlike the Ji — J2 model, this Ising order parameter is 
much less subtle here and the distortions should be more 
prominent. This is exactly what one finds in the phe- 
nomenology of these materials. 

The model presented above gives a simple scenario for 
the observed magnetism and structural transition in the 
pnictide materials. However, it is certainly not a com- 
plete model for the system. At the very least, this phe- 
nomena must coexist with other d electrons, which are 
itinerant and can contribute further to the magnetism. 
Electronic structure studies with local spin-density ap- 
proximation (LSDA) provide further clues to the unusual 
physics. It is found that the local moment and mag- 
netism is highly sensitive to the assumed spin configura- 
tion. Not all choice of LSDA patterns lead to well formed 
moments. The (tt, 0) pattern leads to particularly robust 
magnetism, whereas ferromagnetic ordering leads to al- 
most no moment at all. 

We can view such a calculation as providing further in- 
formation on orbital occupations, when a spin-density is 
imposed on the system. Fixing the (tt, 0)-state in LSDA 
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pattern drives the system to the orbital ordered mag- 
netic state discussed above. In contrast, fixing a fer- 
romagnetic spin configuration, not only need not cause 
orbital polarization, it fails to set the occupation of the 
dxz and dy^ orbitals at a level needed for moment forma- 
tion. This suggests that antiferromagnetic spin arrange- 
ment on neighboring sites, and the resulting gain in the 
exchange energy Ji is necessary to stabilize the local mo- 
ment. One way to view the situation is that electrons in 
other, perhaps weakly correlated, itinerant orbitals act 
as a reservoir for these strongly correlated orbitals. The 
gain in the superexchange energy, when neighboring spins 
are anti-parallel, alters their chemical potential in a way 
as to cause them to become Mott insulating, with occu- 
pancy close to unity and well formed moments. This is 
an example of an orbital selective Mott transition, which 
is driven by superexchange interaction between neighbor- 
ing iron atoms. It would be interesting to further study 
this by Dynamical Mean Field Theory (DMFT) and its 
generalizations!^ The magnetic moment formed in these 
orbitals, via the local Hund's coupling, can further polar- 
ize other d-electrons leading to an even enhanced moment 
as found in the electronic structre calculations. 

An intersting aspect of the model is that while the 
system may be highly frustrated above the ordering tem- 
perature, the spins become unfrustrated below the or- 
dering temperature. Yet, the ordered moment can be re- 
duced significantly by quantum fluctuations because we 
are dealing effectively with a quasi-one dimensional spin- 
half system. Indeed, the electronic structure calculations 
show the bands to have much larger dispersion along 
the direction of antiferromagnetic order Ji being much 
larger than all other exchanges means that there will be 
strong one-dimensional fluctuations and that will reduce 
the ordered moment. Any additional moment in the sys- 
tem, coming from the polarization of other d-electrons 
through the Hunds coupling, will be subservient to this 
primary variable. Hence, many factors can contribute to 
the full moment and indeed these systems show a lot of 
variation from one material to another. 

The most direct way to confirm such a scenario would 
be to look for orbital order through X-ray scattering. 
Since the pioneering studies of Kugel and Khomskii^^ 
the subject of orbital order has come a long way and is 



known to determine key aspects of low energy physics of 
many materials. The electronic structure studies of these 
materials by Han et al imply that even without any lat- 
tice distortion, the tetragonal symmetry must be sponta- 
neously broken. Orbital order is a natural and perhaps 
only candidate for such a broken symmetry, where the ex- 
change constants themselves will be different along the X 
and Y axes. However, almost certainly, the orbital order 
is more subtle in these systems than the simplified model 
suggests. The occupation of d^z orbitals may exceed the 
occupation of dyz orbitals by a small amount. Such a 
small symmetry breaking may be sufficient to cause the 
exchange interactions to be antiferromagnetic in one di- 
rection and ferromagnetic in the other. Nevertheless, our 
model serves as a good zeroth order starting point for un- 
derstanding this behavior. It would be useful to look for 
such occupation numbers in more detail in the electronic 
structure calculations. 

The orbital degrees of freedom should couple strongly 
to the lattice. Hence, they should have long time scales 
associated with them and may act like quenched impuri- 
ties and give rise to glassy behavior. It would be interest- 
ing to look for such effects in experiments. Furthermore, 
the role of vacancies would be very different, in such a sys- 
tem, from usual local moment antifcrromagncts. Rather 
than form moments around vacancies, one is more likely 
to destroy the Mott behavior altogether 

To get a complete picture for this material, one needs 
to consider a multi-band model with all the low energy 
electrons. However, this analysis shows that in some re- 
spects this system may be much closer to the cuprates 
than previously thought and hiding inside the multi- 
electron weakly correlated electronic system maybe a 
spin-half strongly correlated component which drives the 
magnetic phenomena. Whether this also plays a role in 
superconductivity remains to be seen. 
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FIG. 1: Linear spin-wave spectra of the model through the 
Neel point (tt, tt). Case I refers to an antiferromagnetic cou- 
pling along X and a ferromagnetic coupling along y, where as 
case II refers to equal antiferromagnetic coupling along both 
axes. In the latter case, order by disorder phenomena opens 
a gap at {qy = tt) but it leaves the point a mimima. 



